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ABSTRACT 



We review theoretical ideas, problems and implications of neutrino masses and 
mixing angles. We give a general discussion of schemes with three light neutrinos. 
^ ' Several specific examples are analyzed in some detail, particularly those that can 

I/"*) , be embedded into grand unified theories. 

\Q ' 1. Introduction 
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There is by now con vincing evidence, from the experimental study of atmospheric 
and solar neutrinos HE?, for the existence of at least two distinct frequencies of neu- 
trino oscillations. This in turn implies non-vanishing neutrino masses and a mixing 
o ■ matrix, in analogy with the quark sector and the CKM matrix. So apriori the study 

of masses and mixings in the lepton sector should be considered at least as important 
as that in the quark sector. But actually there are a number of features that make 
neutrinos especially interesting. In fact the smallness of neutrino masses is probably 
related to the fact that u's are completely neutral (i.e. they carry no charge which is 
exactly conserved) and are Majorana particles with masses inversely proportional to 
the large scale where lepton number (L) conservation is violated. Majorana masses 
can arise from the see-saw mechanism , in which case there is some relation with the 
Dirac masses, or from higher- dimensional non-renormalizable operators which come 
from a different sector of the lagrangian density than any other fermion mass terms. 
The relation with L non-conservation and the fact that the observed neutrino oscilla- 
tion frequencies are well compatible with a large scale for L non-conservation, points 
to a tantalizing connection with Grand Unified Theories (GUT's). So neutrino masses 
and mixings can represent a probe into the physics at GUT energy scales and offer 
a different perspective on the problem of flavour and the origin of fermion masses. 
There are also direct connections with important issues in astrophysics and cosmology 
as for example baryogenesis through leptogenesis ^ and the possibly non-negligible 
contribution of neutrinos to hot dark matter in the Universe. 



Recently there have been new important experimental results that have consid- 
erably improved our knowledge. The SNO experiment has confirmed that the solar 
neutrino deficit is due to neutrino oscillations and not to a flaw in our modeling of 
the sun: the total neutrino flux is in agreement with the solar model but only about 
one third arrives on earth as u e while the remaining part consists of other kinds of 
active neutrinos, presumably and v T . The allowed amount of sterile neutrinos is 
strongly constrained. The KamLAND experiment has established that FjT from reac- 
tors show oscillations over an average distance of about 180 Km which are perfectly 
compatible with the frequency and mixing angle corresponding to one of the solu- 
tions of the solar neutrino problem (the MSW Large Angle (LA) solution). Thus 
the results from solar neutrinos have been reproduced and improved by a terrestrial 
experiment. Also the coincidence of the frequency for neutrinos from the sun and for 
antineutrinos from reactors is consistent with the validity of CPT invariance, which 
has been questioned because of the puzzling LSND claim of a signal that could in- 
dicate a third distinct oscillation frequency (hence implying either more than three 
light neutrinos or CPT violation). Finally the first results from WMAP on the cosmic 
radiation background and the related determination of cosmological parameters, in 
combination with other measurements, lead to an upper limit on the cosmological 
neutrino density Q u ^ 0.015. This is a very important result that indicates that 
neutrinos are not a major component of the dark matter in the Universe. For three 
degenerate neutrinos the WMAP limit implies an upper bound on the common mass 
given by m u ^ 0.23 eV. Given the priors that are assumed for this determination, 
i.e. a definite cosmological model, a 2-digit value for the bound is not to be taken 
too seriously. Still the quoted value is about an order of magnitude smaller than the 
bound from tritium beta decay and of the same order of the upper bound on the 
Majorana mass that fixes the rate of neutrinoless double beta decay. 

In spite of this progress there are many alternative models of neutrino masses. 
This variety is mostly due to the considerable experimental ambiguities that still 
exist. One first missing input is the absolute scale of neutrino masses: neutrino 
oscillations only determine mass squared differences. For atmospheric neutrinos 
^■ m ltm ~ 2.5 10~ 3 eV 2 while for solar neutrinos Am 2 sol ~ 7 10~ 5 eV 2 . An- 
other key missing quantity is the value of the third mixing angle si 3 on which onlva 
bound is known, S13 ^ 0.22. Then it is essential to know whether the LSND signal'"', 
which has not been confirmed by KARMEN® 1 and is currently being double-checked 
by MiniBoone ^ , will be confirmed or will be excluded. If LSND is right we probably 
need at least four light neutrinos; if not we can do with only the three known ones. 

Here we will briefly summarize the main categories of neutrino mass models, dis- 
cuss their respective advantages and difficulties and give a number of examples. We 
illustrate how forthcoming experiments can discriminate among the various alterna- 
tives. We will devote a special attention to a comprehensive discussion in a GUT 
framework of neutrino masses together with all other fermion masses. This is for ex- 



ample possible in models based on SU(5)x U(l)p or on SO(10) (we always consider 
SUSY GUT's)E0. 



2. Basic Formulae and Data for Three-Neutrino Mixing 

We assume in the following that the LSND signal^, will not be confirmed so that 
there are only two distinct neutrino oscillation frequencies, the atmospheric and the 
solar frequencies. These two can be reproduced with the known three light neutrino 
species (for more than three neutrinos see, for example, ref. ^^). 

Neutrino oscillations are due to a misalignment between the flavour basis, v' = 
{y e , i/p, u T ), where v e is the partner of the mass and flavour eigenstate e~ in a left- 
handed (LH) weak isospi n SU(2 ) doublet (similarly for and z/ r )) and the mass 
eigenstates v = z/ 2 , z/ 3 ) IHUST 

v' = Uu , (1) 

where U is the unitary 3 by 3 mixing matrix. Given the definition of U and the 
transformation properties of the effective light neutrino mass matrix m v : 

v' T m v v' = v T U T m v Uv (2) 
U T m u U = Diag (mi, m 2 ,m 3 ) = m diag , 

we obtain the general form of m v (i.e. of the light v mass matrix in the basis where 
the charged lepton mass is a diagonal matrix): 

m v = Um diag U T . (3) 

The matrix U can be parameterized in terms of three mixing angles #12, #23 and #13 
(0 < 9ij < 7r/2) and one phase (p (0 < ip < 2ti) exactly as for the quark mixing 
matrix Vckm- The following definition of mixing angles can be adopted: 



U 




( c 13 s 13 e^ 

1 I [ -s 12 c 12 I (4) 
\-s 13 e~^ 




where sy = sin^-, = cos%. In addition, if v are Majorana particles, we have the 
relative phases among the Majorana masses m 1; m 2 and ma. If we choose m 3 real 
and positive, these phases are carried by m 12 = Im^le'^ 1 ' 2 Thus, in general, 9 
parameters are added to the SM when non-vanishing neutrino masses are included: 
3 eigenvalues, 3 mixing angles and 3 CP violating phases. 

In our notation the two frequencies, Am'j/AE (I = sun,atm), are parametrized 
in terms of the v mass eigenvalues by 

A ™Ln = l Am ? 2 | 5 Am ltm = l Am 2 3 l • ( 5 ) 



Table 1: Square mass differences and mixing angles . 
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where Am 2 2 = |m 2 | 2 — Im^ 2 > and Am 2 3 = m\ — |m 2 | 2 . The numbering 1,2,3 
corresponds to our definition of the frequencies and in principle may not coincide 
with the ordering from the lightest to the heaviest state. From ex peri ment, see table 
[TJ we know that si 3 is small, according to CHOOZ, s Vi < 0.22 (3a) Atmospheric 

neutrino oscillations mainly depend on (Am 2 tm , #23, #13), while solar oscillations are 
controlled by (Am 2 oi , #i 2 , 13 ). Therefore, in the ideal limit of exactly vanishing S13, 
the solar and atmospheric oscillations decouple and depend on two separate sets 
of two-flavour parameters. For atmospheric neutrinos we have C23 ~ S23 ~ 1/ \/2, 
corresponding to nearly maximal mixing. Oscillations of muon neutrinos into tau 
neutrinos are favoured over oscillations into sterile neutrinos (v s ). The conversion 
probability and the zenith angular distribution of high-energy muon neutrinos are 
sensitive to matter effects, which distinguish v T from u s . Moreover, for conversion 
of into pure u s , neutral current events would become up/down asymmetric. In 
both cases data strongly disfavour the pure sterile case. Oscillations into v T are also 
indirectly supported by a SK data sample that can be interpreted in terms of enriched 
r-like charged-current events. The sterile component of the neutrino participating in 
atmospheric oscillations should square to less than 0.25 at 90% CL. Disappearance 



of laboratory-produced muon neutrinos has also been confirmed within expectations 
by the K2K experiment. 

The only surviving solution to the solar neutrino problem after KamLAND results 
is LA MSW, with Am 2 sol « 7 • 1CT 5 eV 2 and tan 2 12 ~ 0.4. At the moment data are 
comp atible with two small 'islands', almost separate along Am 2 oi (see table 1, from 
I I '), with a slight preference for the low Am 2 o( region. Before KamLAND the 
interpretation of solar neutrino data in terms of oscillations required the knowledge 
of the Boron neutrino flux, For instance, charged and neutral current data from 
SNO are sensitive, respectively, to /s(-P ee ) and to fB(J2 a Pea) ( a = e, u, r), where (P e f) 
denotes the appropriately averaged conversion probability from v e to i/f. KamLAND 

provides a direct measurement of (P ee )- Beyond the impact on the oscillation 
parameters and a check that the solar standard model works well (fg = (1.00 ± 
0.06) x 5.05 • 10 6 cm~ 2 sec -1 ) , the comparison among these experiments shows that 
the conversion of Boron solar neutrinos into sterile neutrinos is compatible with zero. 
Therefore, the LSND indication for a third oscillation frequency associated to one 
or more sterile neutrinos is not supported by any other experiment, at the moment. 
Now, after KamLAND, also the possibility that such a frequency originates from a 
CPT violating neutrino spectrum has no independent support. Data from solar 
neutrino experiments and KamLAND, involving, respectively, electron neutrinos and 
electron antineutrinos, are compatible with a CPT invariant spectrum. 

If we take maximal S23 and keep only linear terms in u = si^e tip , from experiment 
we find the following structure of the Up (/ = e,/i,r, i — 1, 2, 3) mixing matrix, apart 
from sign convention redefinitions: 

/ C12 _ S12 _ u \ 

U fl = \ -{s 12 + c 12 u*)/V2 (c 12 -s 12 u*)/V2 1/V2 , (6) 
V (sia - c 12 u*)/V2 -(c 12 + s 12 u*)/V2 1/^2 J 

where 9\ 2 is close to n/6. Given the observed frequencies and our notation in eq. (jSJ, 
there are three possible patterns of mass eigenvalues: 



Degenerate 
Inverted hierarchy 
Normal hierarchy 



I m-i I ~ \m 2 \ ~ |m 3 | ^> {rrii — mj\ 
\ m i\ ~ \ m 2\ 3> 1 77l3| 

\m 3 \ > |m 2 ,i| (7) 



Models based on all these patterns have been proposed and studied and all are in 
fact viable at present. In the following we will first discuss neutrino masses in general 
and, in particular, Majorana neutrinos. Then we recall the existing constraints on 
the absolute scale of neutrino masses. We then discuss the importance of neutrinoless 
double beta decay that, if observed, would confirm the Majorana nature of neutrinos. 
Also the knowledge of the rate of this process could discriminate among the possi- 
ble patterns of neutrino masses in (J7J). The possible importance of heavy Majorana 



neutrinos for the explanation of baryogenesis through leptogenesis in the early Uni- 
verse will be briefly discussed. We finally review the phenomenology of neutrino mass 
models based on the three spectral patterns in (JZJ) and the respective advantages and 
problems. 



3. Neutrino Masses and Lepton Number Violation 

Neutrino oscillations imply neutrino masses which in turn demand either the ex- 
istence of right-handed (RH) neutrinos (Dirac masses) or lepton number L violation 
(Majorana masses) or both. Given that neutrino masses are certainly extremely 
small, it is really difficult from the theory point of view to avoid the conclusion that 
L conservation must be violated. In fact, in terms of lepton number violation the 
smallness of neutrino masses can be explained as inversely proportional to the very 
large scale where L is violated, of order M GUT or even M Pl . 

Once we accept L non- conservation we gain an elegant explanation for the small- 
ness of neutrino masses. If L is not conserved, even in the absence of heavy RH 
neutrinos, Majorana masses for neutrinos can be generated by dimension five opera- 
tors of the form 

o 5 = iMMElk + h . c . j (8) 

with H being the ordinary Higgs doublet, U the SU(2) lepton doublets, A a matrix in 
flavour space and A a large scale of mass, of order Mqut or Mpi. Neutrino masses 
generated by O5 are of the order m u m v 2 /A for Ay ~ 0(1), where v ~ 0(100 GeV) 
is the vacuum expectation value of the ordinary Higgs. 

We consider that the existence of RH neutrinos u c is quite plausible because all 
GUT groups larger than SU(5) require them. In particular the fact that u c completes 
the representation 16 of SO(10): 16=5+10+1, so that all fermions of each family are 
contained in a single representation of the unifying group, is too impressive not to 
be significant. At least as a classification group SO(10) must be of some relevance. 
Thus in the following we assume that there are both v c and L non- conservation. With 
these assumptions the see-saw mechanism^ is possible. Also to fix notations we recall 
that in its simplest form it arises as follows. Consider the SU(3) x SU(2) x U(l) 
invariant Lagrangian giving rise to Dirac and v c Majorana masses (for the time being 
we consider the v (versus v c ) Majorana mass terms as comparatively negligible): 

£ = -v cT y u (Hl) + V T Mz/ c + h.c. (9) 

The Dirac mass matrix m^E y v v / y/2, originating from electroweak symmetry break- 
ing, is, in general, non-hermitian and non- symmetric, while the Majorana mass matrix 
M is symmetric, M = M T . We expect the eigenvalues of M to be of order Mqut 
or more because v c Majorana masses are SU(3)x SU(2)x U(l) invariant, hence un- 
protected and naturally of the order of the cutoff of the low-energy theory. Since all 



v c are very heavy we can integrate them away. For this purpose we write down the 
equations of motion for v c in the static limit, i.e. neglecting their kinetic terms: 



- — = y u (Hl) - Mu c = . (10) 

From this, by solving for z/ c , we obtain: 

u c = M~ l y u {Hl) . (11) 

We now replace in the lagrangian, eq. ©, this expression for v c and we get the 
operator O5 of eq. (JEJ) with 

^ = -vlM-\ v , (12) 

and the resulting neutrino mass matrix reads: 

m v — Tr!F D M~ x mD . (13) 

This is the well known see-saw mechanism result ^ : the light neutrino masses are 
quadratic in the Dirac masses and inversely proportional to the large Majorana mass. 
If some v c are massless or light they would not be integrated away but simply added 
to the light neutrinos. Notice that the above results hold true for any number n of 
heavy neutral fermions R coupled to the 3 known neutrinos. In this more general case 
M is an n by n symmetric matrix and the coupling between heavy and light fields is 
described by the rectangular n by 3 matrix mo- Note that for m v m J Aml tm ~ 0.05 
eV and m v « m 2 D /M with mo ~ v « 200 GeV we find M w 10 15 GeV which indeed 
is an impressive indication for Mqjjt- 

If additional non-renormalizable contributions to O5, eq. |Hl are comparatively 
non- negligible, they should simply be added. After elimination of the heavy right- 
handed fields, at the level of the effective low-energy theory, the two types of terms 
are equivalent. In particular they have identical transformation properties under a 
chiral change of basis in flavour space. The difference is, however, that in the see-saw 
mechanism, the Dirac matrix mo is presumably related to ordinary fermion masses 
because they are both generated by the Higgs mechanism and both must obey GUT- 
induced constraints. Thus if we assume the see-saw mechanism more constraints are 
implied. 

4. Importance of Neutrinoless Double Beta Decay 

Oscillation experiments do not provide information about the absolute neutrino 
spectrum and cannot distinguish between pure Dirac and Majorana neutrinos. From 
the endpoint of tritium beta decay spectrum we have an absolute upper limit of 2.2 



eV (at 95% C.L.) on the mass of electron antineutrino , which, combined with the 
observed oscillation frequencies under the assumption of three CPT-invariant light 
neutrinos, represents also an upper bound on the masses of the other active neutri- 
nos. Complementary information on the sum of neutrino masses is also provided by 
the galaxy power spectrum combined with measurements of the cosmic microwave 
background anisotropies. According to the recent analysis of the WMAP collabo- 
ration^^, J2i \ m i\ < 0-69 eV (at 95% C.L.). More conservative analyses -23' give 
J2i \ m i\ < 101 eV, still much more restrictive than the laboratory bound. 

The discovery of QvfiP decay would be very important because it would establish 
lepton number violation and the Majorana nature of z/s, and provide direct informa- 
tion on the absolute scale of neutrino masses. As already mentioned the present limit 
from 0vj3j3 is \m ee \ < 0.2 eV or to be more conservative \m ee \ < 0.3 — 0.5 eV^^. 
Note, however, that the WMAP limit implies for 3 degenerate z/s \m\ ^ 0.23 eV and 
this poses a direct constraint on m ee . 

It is interesting to see what is the level at which a si gnal c an be expected or at 
least not excluded in the different classes of models in (J7J) ESESJ The quantity which 
is bound by experiments is the 11 entry of the v mass matrix, which in general, from 
eqs. (J2J) and (HJ), is given by : 

\m ee \ = |(1- s 2 13 ) {m x cl 2 + m 2 s 2 12 ) + m 3 e 2i ^s 2 13 \ , (14) 

For 3-neutrino models with degenerate, inverse hierarchy or normal hierarchy mass 
patterns, starting from this general formula it is simple to derive the following bounds. 

a) Degenerate case. If \m\ is the common mass and we take S\ 3 = 0, which is a safe 
approximation in this case, because \m 3 \ cannot compensate for the smallness 
of s\s, we have m ee ~ |m|(cf 2 ±s^ 2 ). Here the phase ambiguity has been reduced 
to a sign ambiguity which is sufficient for deriving bounds. So, depending on the 
sign we have m ee = \m\ or m ee = \m\cos26\2- We conclude that in this case m ee 
could be as large the present experimental limit but should be at least of order 
0(\J Am 2 trn ) ~ O(10~ 2 eV) unless the solar angle is practically maximal, in 
which case the minus sign option can be arbitrarily small. But the experimental 
l-cr range of the solar angle does not favour a cancellation by more than a factor 
of 3. 

b) Inverse hierarchy case. In this case the same approximate formula m ee = 
\m\(c 2 2 ±s 2 2 ) holds because m 3 is small and si 3 can be neglected. The difference 
is that here we know that \m\ « J Am 2 tm so that \m ee \ < \J Am 2 tm ~ 0.05 
eV. 

c) Normal hierarchy case. Here we cannot in general neglect the m 3 term. However 
in this case \m ee \ ~ | ^ Am 2 sun s\ 2 ± yj Am 2 tm s 2 3 \ and we have the bound 
\m ee \ < a few 10 -3 eV. 



Recently evidence for 0vf3f3 was claimed in ref. at the 2-3<r level corresponding 
to |m ee | ~ 0.39 eV. If confirmed this would rule out cases b) and c) and point to 
case a) or to models with more than 3 neutrinos. 

5. Baryogenesis via Leptogenesis from Heavy v c Decay 

In the Universe we observe an apparent excess of baryons over antibaryons. It is 
appealing that one can explain the observed baryon asymmetry by dynamical evo- 
lution (baryogenesis) starting from an initial state of the Universe with zero baryon 
number. For baryogenesis one needs the three famous Sakharov conditions: B viola- 
tion, CP violation and no thermal equilibrium. In the history of the Universe these 
necessary requirements can have occurred at different epochs. Note however that 
the asymmetry generated by one epoch could be erased at following epochs if not 
protected by some dynamical reason. In principle these conditions could be verified 
in the SM at the electroweak phase transition. B is violated by instantons when kT 
is of the order of the weak scale (but B-L is conserved), CP is violated by the CKM 
phase and sufficiently marked out-of- equilibrium conditions could be realized during 
the electroweak phase transition. So the conditions for baryogenesis at the weak scale 
in the SM superficially appear to be present. However, a more quantitative analysis 

shows that baryogenesis is not possible in the SM because there is not enough 
CP violation and the phase transition is not sufficiently strong first order, unless 
mfl < 80 GeV, which is by now completely excluded by LEP. In SUSY extensions 
of the SM, in particular in the MSSM, there are additional sources of CP violation 
and the bound on rrin is modified by a sufficient amount by the presence of scalars 
with large couplings to the Higgs sector, typically the s-top. What is required is that 
ruh ~ 80 — 110 GeV, a s-top not heavier than the top quark and, preferentially, a 
small tan/3. But also this possibility has by now become at best marginal with the 
results from LEP2. 

If baryogenesis at the weak scale is excluded by the data it can occur at or just 
below the GUT scale, after inflation. But only that part with |B — L| > would 
survive and not be erased at the weak scale by instanton effects. Thus baryogenesis 
at kT ~ 10 10 — 10 15 GeV needs B-L violation at some stage like for m u if neutrinos 
are Majorana particles. The two effects could be related if baryogenesis arises from 
leptogenesis then converted into baryogenesis by instantons ™. Recent results on 
neutrino masses are compatible with this elegant possibility Thus the case of 
baryogenesis through leptogenesis has been boosted by the recent results on neutrinos 

In leptogenesis the departure from equilibrium is determined by the deviation 
from the average number density induced by the decay of the heavy neutrinos. The 
Yukawa interactions of the heavy Majorana neutrinos v c lead to the decays v c — ► IH 
(with / a lepton) and v c — > IH with CP violation. The violation of L conservation 



arises from the AL = 2 terms that produce the Majorana mass terms. The rates of 
the various interaction processes involved are temperature dependent with different 
powers of T, so that the equilibrium densities and the temperatures of decoupling 
from equilibrium during the Universe expansion are different for different particles 
and interactions. The rates T& L (T) of AL = 2 processes depend also on the neu- 
trino masses and mixings, so that the observed values of the baryon asymmetry are 
related to neutrino processes. A recent quantitative analysis in - ■ concludes that 
a successful baryogenesis requires neutrino masses in the range ^ 0.1 eV. This 
range is remarkably consistent with the results on neutrino oscillations. Beyond ther- 
mal production after reheating, RH neutrinos can also originate from large inflaton 
oscillations during the preheating stage 



6. Degenerate Neutrinos 

For degenerate neutrinos the average m 2 is much larger than the splittings. At first 
sight the degenerate case is the most appealing: the observation of nearly maximal 
atmospheric neutrino mixing and the more recent result that also the solar mixing 
is large suggests that all v masses are nearly degenerate. We shall see that this 
possibility has become less attractive with the recent new experimental information. 

It is clear that in the degenerate case the most likely origin of v masses is from 
some dimension 5 operators (Hl)J Xij(Hl)j / A not related to the see-saw mechanism 
m v = mpM^rriD- In fact we expect the v Dirac mass mo not to be degenerate like for 
all other fermions and a conspiracy to reinstate a nearly perfect degeneracy between 
mo and M, which arise from completely different physics, looks very unplausible (see, 
however, Thus in degenerate models, in general, there is no direct relation with 
Dirac masses of quarks and leptons and the possibility of a simultaneous description 
of all fermion masses within a grand unified theory is more remote . 

The degeneracy of neutrinos should be guaranteed by some slightly broken sym- 
metry. Models based on discrete or continuous symmetries have been proposed. For 
example in the models of ref. the symmetry is SO (3): in the unbroken limit 
neutrinos are degenerate and charged leptons are massless. When the symmetry is 
broken the charged lepton masses are much larger than neutrino splittings because the 
former are first order while the latter are second order in the electroweak symmetry 
breaking. 

The upper limit on the common value \m\ becomes particularly stringent if one 
adopts the cosmological WMAP bound \m\ ^ 0.23 eV. The more direct laboratory 
limit from tritium beta decay is \m\ ^ 2.2 eV. In past years degenerate models 
with v masses as large as \m\ ~ 1 — 2 eV were considered with the perspective 
of a large fraction of hot dark matter in the universe. In this case, however, the 
existing limit on the absence of 0z//3/3 (|m ee | < 0.2 eV or to be more conservative 
|m ee | < 0.3 — 0.5 eV) implies approximate double maximal mixing (bimixing) for 



solar and atmospheric neutrinos. As discussed in sect. 4 , for |m| » m ee , one needs 
mi ~ —^2 and, to a good accuracy, c 2 2 ~ s 2 2 , in order to satisfy the bound on m ee . 
This is exemplified by the following texture 

/ -l/y/2 l/y/2 \ 

m u = m\ -l/y/2 (l + r])/2 (1 + ^/2 , (15) 
V l/y/2 (1 + 77)/2 (l + v)/2j 

where r|<l, corresponding to an exact bimaximal mixing, s 13 = and the eigenval- 
ues are mi = m, m 2 = — m and m 3 = (1+77)771. This texture has been proposed in the 
context of a spontaneously broken SO (3) flavor symmetry and it has been stu died to 
analyze the stability of the degenerate spectrum against radiative corrections l ^ | 38 ) 
A more realistic mass matrix can be obtained by adding small perturbations to m u 
in eq. (fTHJ) : 

/ 5 -l/y/2 (l-e)/V2 \ 

m u = m\ -l/y/2 (1 + V )/2 (l + r)-e)/2 , (16) 

V(l-e)/V2 (l + tj-e)/2 (1 + 77 - 2e)/2 7 

where e parametrizes the leading flavor- dependent radiative corrections (mainly in- 
duced by the r Yukawa coupling) and 5 controls m ee . Consider first the case 5 ^ e. 
To first approximation 12 remains maximal. We get Am? sun w m 2 e 2 /r] and 

/Am 2 \ 1/2 /Am 2 Am 2 \ 1/2 

/ ^ , m ee < m gtm f ww . (17) 



713 ~ I Am 2 



aim , 



If we instead assume 5 > e, we find Am 2 nn 2m 2 5, 23 ~ 7r/4, sin 2 26*i 2 ~ 1 — 5 2 /4. 
Also in this case the solar mixing angle remains close to 7r/4. We get: 



Am 2 

#13~0 , /,/, ' M "' 



2m 

too small for detection if the average neutrino mass m is around the eV scale. We see 
that with increasing \m\ more and more fine tuning is needed to reproduce the LA 
solution values of Am 2 un and #12- In conclusion, even without invoking the WMAP 
limit, large v masses, |m| ~ 1 — 2 eV, are disfavoured by the limit on Qu(3(3 decay 
and by the emerging of the LA solution with the solar angle definitely not maximal. 
Also, we have seen in sect. 5 that the attractive mechanism of baryogenesis through 
leptogenesis appears to disfavour \m\ ^ 0.1 eV. All together, after WMAP and 
KamLAND, among degenerate models those with |m| ^ 0.23 eV are favoured by 
converging evidence from different points of view. 

For \m\ not larger than the Qv/3/3 bound one does not need a cancellation in m ee 
and mi and m 2 can be approximately equal in magnitude and phase. For example, 



in the limit S13 = 0, the matrix 



m. 



m 



( 1 










(19) 







corresponds to maximal #23 (pseudo Dirac 23 sub- matrix) with Diag[m v ] — m(l, 1, — 1) 
The angle 9 12 is unstable and a small perturbation can give any value to it. Note, 
however, that in this case the non vanishing matrix elements must be of equal abso- 
lute value and not just of order 1. So either this is guaranteed by a symmetry (as, 



for example, in ref. 



35, 



or the model is unnaturally fine-tuned. 



As a different example (also with no cancellation between mi and m%), a model, 
which is simple to describe but difficult to derive in a natural way, is one where 
up quarks, down quarks and charged leptons have "democratic" mass matrices, with 
all entries equal (in first approximation): 



rrif = rrif 




+ 5rrif 



(20) 



where rhf (f = u,d,e) are three overall mass parameters and 5mf denote small 
perturbations. If we neglect 5m f, the eigenvalues of rrif are given by (0, 0, 3 rhf). The 
mass matrix rrif is diagonalized by a unitary matrix Uf which is in part determined 
by the small term 8m f. If 5m u « 5m<i, the CKM matrix, given by Vckm = U}jJd, is 
nearly diagonal, due to a compensation between the large mixings contained in U u and 
lid- When the small terms 5mf are diagonal and of the form 5mf = Diag(— 6f, 6f, 5f) 
the matrices U f are approximately given by (note the analogy with the quark model 
eigenvalues n°, rj and rj'): 



1/V2 -1/V2 
v) & I 1/V6 i/VE -2/V6 

l/y/3 1/V3 l/y/3 



(21) 



At the same time, the lightest quarks and charged leptons acquire a non-vanishing 
mass. The leading part of the mass matrix in eq. (J2*U|) is invariant under a discrete 
S3L x 1S3.R permutation symmetry. The same requirement leads to the general neutrino 
mass matrix: 







f 1 
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f 1 
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m v = m 
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v0 









V 1 
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5m, 



(22) 



where 5m v is a small symmetry breaking term and the two independent invariants 
are allowed by the Majorana nature of the light neutrinos. If r vanishes the neutrinos 
are almost degenerate. In the presence of 5m v the permutation symmetry is broken 



and the degeneracy is removed. If, for example, we choose 5m u = Diag(0, e, rj) with 
e < ») < 1 and r e, the solar and the atmospheric oscillation frequencies are 
determined by e and 77, respectively. The mixing angles are almost entirely due to 
the charged lepton sector. The problems with this class of models is that the solar 
angle should be close to maximal (typically more than the atmospheric angle) and 
that the neutrino spectrum and mixing angles are not determined by the symmetric 
limit, but only by a specific choice of the parameter r and of the perturbations that 
cannot be easily justified on theoretical grounds. 

In conclusion, the parameter space for degenerate models has recently become 
smaller because of the indications from WMAP and from leptogenesis that tend to 
lower the maximum common mass allowed for light neutrinos. It is also rather difficult 
to reproduce the observed pattern of frequencies and mixing angles, in particular 
two large and one small mixing angle, with the solar angle large but not maximal. 
Degenerate models that fit can only arise from a very special dynamics or a non 
abelian flavour symmetry with suitable breakings. 

6.1. Anarchy 

Anarchical models can be considered as particular cases of degenerate models 
with m 2 ~ Am 2 tm . In this class of models mass degeneracy is replaced by the principle 
that all mass matrices are structure-less in the neutrino (or even in the whole lepton) 
sector. For the MSW LA solution the ratio of the solar and atmospheric frequencies 
is not so small: r = (Aml un ) LA / Am^ tm ~ 1/40 and two out of three mixing angles 
are large. The key observation is that the see-saw mechanism tends to enhance the 
ratio of eigenvalues: it is quadratic in mo so that a hierarchy factor /inmo becomes 
P in m v and the presence of the Majorana matrix M results in a further widening 
of the distribution. Another squaring takes place in going from the masses to the 
oscillation frequencies which are quadratic. As a result a random generation of the 
mo and M matrix elements leads to a distribution of r that peaks around 0.1. At 
the same time the distribution of sin 2 9y is rather flat for all three mixing angles. 
Clearly the smallness of #13 is problematic for anarchy. This can be turned into the 
prediction that in anarchical models #13 must be near the present bound (after all the 
value 0.2 for sin # 13 is not that smaller than the maximal value 0.701). In conclusion, 
if #13 is near the present bound, one can argue that the neutrino masses and mixings, 
interpreted by the see-saw mechanism, can just arise from structure-less underlying 
Dirac and Majorana matrices. 

7. Inverted Hierarchy 

The inverted hierarchy configuration \m\\ ~ \m2\ ^> \rris\ consists of two levels 
mi and 7712 with small splitting Am 2 2 = Am 2 sun and a common mass given by 



\ m h2\ 



\ Am atm\ 



3 • 10~ 3 eV 2 . One particularly interesting example of this sort^^, 
which leads to double maximal mixing, is obtained with the phase choice mi = — m 2 
so that, approximately: 



m diag = Diag(v / 2m, -y/2m, 0) . 
The effective light neutrino mass matrix 



m,. 



Um diag U 



(23) 



(24) 



which corresponds to the mixing matrix of double maximal mixing cyi — syi — 1/ \/2 
and s\3 = u = in eq. (JHJ). 



E7 



/ l/y/2 l/y/2 

-1/2 1/2 l/y/2 
\ 1/2 -1/2 1^/2, 



(25) 



is given by: 

/ -1 l\ 
m„ = m -1 . (26) 
V 1 0/ 

The structure of m u can be reproduced by imposing a U(l) flavour symmetry with 
charge L e — — L r starting either from (Hl)J\ij(Hl)j/A or from RH neutrinos 
via the see-saw mechanism. The 1 — 2 degeneracy remains stable under radiative 
corrections. 

The leading texture in (|2*6^) can be perturbed by adding small terms: 



77V 




(27) 



where 5 and rj are small (<^C 1), real parameters defined up to coefficients of order 
1 that can differ in the various matrix elements. One could also make the absolute 
values of the 12, 13 terms different by terms of order s^S. The perturbations leave 
/S.m 2 atm and #23 unchanged, in first approximation. We obtain tan 2 9 12 ~ 1 + S + rj and 
Am 2 un / /S.m 2 atm pa T] + 8, where coefficients of order one have been neglected. Moreover 



#13 ~ rj. If rj 3> 5, we have 



'13 



Am lan 
Am atm 



m ee < 



Am2 sun I A/;r> 



(2? 



atm , 



In the other case, 77 <C 5 we obtain: 



»<*ee 2 V sun 



yAm 2 tm 



(29) 



There is a well-known difficulty of this scenario to fit the MSW LA solution | 41 | 42 | 43 | 
Indeed, barring cancellation between the perturbations, in order to obtain a Am^ un 
close to the best fit MSW LA value, 77 and 5 should be smaller than about 0.1 and 
this keeps the value of sin 2 29u very close to 1, somewhat in disagreement with global 
fits of solar data^T Even by allowing for a Am 2 nra in the upper range of the MSW 
LA solution, or some fine-tuning between rj and S, we would need large values of the 
perturbations to fit the MSW LA solution values of Am 2 un and 6>i 2 . However, in a 
subclass of models the pattern of parameters needed to bring the solar angle down 
from the maximal value is obtained in a natural way as an effect of the charged lepton 
matrix diagonalisation 

With the phase choice m x = m 2 , i. e. for Diag[m u ] = 771(1, 1,0), in the limit 
S13 = 0, one obtains the matrix 



m u = m 1/2 -1/2 I , (30) 




which corresponds to large #23 with the solar angle unstable to small perturbations. 
In this case fine tuning or a symmetry is necessary to fix the ratios of the matrix 
elements as indicated. 

In conclusion, also for inverse hierarchy some special dynamics or symmetry is 
needed to reproduce the observed features of the data. 



8. Normal Hierarchy 

We now discuss the class of models which we consider the simplest approach to 
neutrino masses and mixings. In particular, in this context one can formulate the most 
constrained framework which allows a comprehensive combined study of all fermion 
masses in GUT's. We start by assuming three widely split i/'s and the existence of 
a RH neutrino for each generation, as required to complete a 16 dimensional repre- 
sentation of SO (10) for each generation. We then assume dominance of the see-saw 
mechanism m u = mF D M~ x rnY). We know that the third-generation eigenvalue of the 
Dirac mass matrices of up and down quarks and of charged leptons is systematically 
the largest one. It is natural (although not necessary) to imagine that this property 
could also be true for the Dirac mass of z/'s: m D U19 ~ Diag(0, 0, muz)- After see-saw 
we expect m v to be even more hierarchical being quadratic in (barring fine-tuned 
compensations between mu and M). Note however tha the amount of hierarchy, 
r = Am^ im /Amj m = m\jm\ is moderate for the MSW LA solution: r ~ 1/40. 

A possible difficulty for hierarchical models is that one is used to expect that 
large splittings correspond to small mixings because normally only close-by states are 
strongly mixed. In a 2 by 2 matrix context the requirement of large splitting and 
large mixings leads to a condition of vanishing determinant and large off-diagonal 



elements. For example the matrix 



x 2 x 



(31) 



has eigenvalues and 1 + x 2 and for x of 0(1) the mixing is large. Thus in the 
limit of neglecting small mass terms of order m 12 the demands of large atmospheric 
neutrino mixing and dominance of m% translate into the condition that the 2 by 2 
sub determinant 23 of the 3 by 3 mixing matrix approximately vanishes. The problem 
is to show that this vanishing can be arranged in a natural way without fine tuning. 
Once near maximal atmospheric neutrino mixing is reproduced the solar neutrino 
mixing can be arranged to be either small or large without difficulty by implementing 
suitable relations among the small mass terms. 

It is not difficult to imagine mechanisms that nat urally lead to the approximate 
vanishing of the 23 sub-determinant. For example | 45 | 46 t assum es that one v c is 
particularly light and coupled to \x and r. In a 2 by 2 simplified context if we have 

»<*(l ?) ■ ^MT o) ■ -»»=(: 5) • < 32 > 

then for a generic mo we find 

m v = m T D M^m D ^ l -(^ h . (33) 

A different possibility that we find attractive is that, in the limit of neglecting terms 
of order my and, in the basis where charged leptons are diagonal, the Dirac m atrix 
mo, defined by v c mov, takes the approximate form, called "lopsided" | 47 | 48 | 49 |i . 



m D oc . (34) 




This matrix has the property that for a generic Majorana matrix M one finds: 

/0 0\ 

m v = m T D M~ x mo oc x 2 x \ . (35) 

\0 x 1 J 

The only condition on M~ l is that the 33 entry is non zero. However, when the 
approximately vanishing matrix elements are replaced by small terms, one must also 
assume that no new 0(1) terms are generated in m v by a compensation between small 
terms in and large terms in M _1 . It is important for the following discussion to 
observe that given by eq. (|34|) under a change of basis transforms as mj) ^ 
V^mnU where V and U rotate the right and left fields respectively. It is easy to 
check that in order to make mo diagonal we need large left mixings (i.e. large off 



diagonal terms in the matrix that rotates LH fields). Thus the question is how to 
reconcile large LH mixings in the leptonic sector with the observed near diagonal 
form of Vckm, the quark mixing matrix. Strictly speaking, since Vckm = U^Ud, 
the individual matrices U u and Ua need not be near diagonal, but Vckm does, while 
the analogue for leptons apparently cannot be near diagonal. However for quarks 
nothing forbids that, in the basis where m u is diagonal, the d quark matrix has large 
non diagonal terms that can be rotated away by a pure RH rotation. We suggest that 
this is so and that in some way RH mixings for quarks correspond to LH mixings for 
leptons. 

In the context of (SUS Y) SU( 5) there is a very attractive hint of how the present 
mechanism can be realized EHED. i n the 5 of SU(5) the d c singlet appears together 
with the lepton doublet (u, e). The (u,d) doublet and e c belong to the 10 and v c to 
the 1 and similarly for the other families. As a consequence, in the simplest model 
with mass terms arising from only Higgs pentaplets, the Dirac matrix of down quarks 
is the transpose of the charged lepton matrix: = {m{) T . Thus, indeed, a large 
mixing for RH down quarks corresponds to a large LH mixing for charged leptons. 
At leading order we may have the lopsided texture: 



In the same simplest approximation with 5 or 5 Higgs, the up quark mass matrix 
is symmetric, so that left and right mixing matrices are equal in this case. Then 
small mixings for up quarks and small LH mixings for down quarks are sufficient 
to guarantee small Vckm mixing angles even for large d quark RH mixings, lepton 
matrices are exactly the transpose of one another cannot be exactly true because of 
the e/d and \ijs mass ratios. It is also known that one remedy to this problem is to 
add some Higgs component in the 45 representation of SU(5) But the symmetry 
under transposition can still be a good guideline if we are only interested in the order 
of magnitude of the matrix entries and not in their exact values. Similarly, the Dirac 
neutrino mass matrix mc is the same as the up quark mass matrix in the very crude 
model where the Higgs pentaplets come from a pure 10 representation of SO(10): 
rriD = m u . For the dominance of the third family eigenvalue as well as a near 
diagonal form could be an order of magnitude remnant of this broken symmetry. 
Thus, neglecting small terms, the neutrino Dirac matrix in the basis where charged 
leptons are diagonal could be directly obtained in the form of eq. (j34jl . 

To get a realistic mass matrix, we allow for deviations from the symmetric limit 
of ()35|) with x = 1. For instance, we can consider those models where the neutrino 
mass matrix elements are dominated, via the see-saw mechanism, by the exchange of 
two right-handed neutrinos z ^'. Since the exchange of a single RH neutrino gives a 
successful zeroth order texture, we are encouraged to continue along this line. Thus, 




/0 
1 
\0 1 




(36) 



we add a sub-dominant contribution of a second RH neutrino, assuming that the third 
one gives a negligible contribution to the neutrino mass matrix, because it has much 
smaller Yukawa couplings or is much heavier than the first two. The Lagrangian that 
describes this plausible subset of see-saw models, written in the mass eigenstate basis 
of RH neutrinos and charged leptons, is 



£ = y i u c Hl i + y ym i + ^v c2 + ^-v c ' 2 , (37) 



leading to 



x ViVj . ViVj foe 



M M- 

where i,j = {e, /i, r}. In particular, if y e <C y^ pa y T and y' ~ y' T , we obtain: 



m,, = m 



( 6 

e 1+ry 1+v | , (39) 




where coefficients of order one multiplying the small quantities 5, e and r\ have been 
omitted. The 23 sub determinant is generically of order rj. The mass matrix in (|39p 
does not describe the most general perturbation of the zeroth order texture We 
have implicitly assumed a symmetry between and v T which is preserved by the 
perturbations, at least at the level of the order of magnitudes. The perturbed texture 
(}3*9*j) can also arise when the zeroes of the lopsided Dirac matrix in (}3"4*j) are replaced 
by small quantities. It is possible to construct models along this line based on a 
spontaneously broken U(l)p flavor symmetry, where 5, e and rj are given by positive 
powers of one or more symmetry breaking parameters. Moreover, by playing with 
the U ( 1 ) I p charges , we can adjust, to certain extent, the relative hierarchy between r/, 
e and S | 45 | 46 | 48 | 49 | 50 | ^3 )^ ag we w jjj gee j n sec ^i on g_ The texture (}3*9*j) can also be 

generated in SUSY models with i?-parity violation 

Let us come back to the mass matrix m v of eq. (|39p. After a first rotation by an 
angle 9 2 3 close to 7r/4 and a second rotation with 6> 13 pa e, we get 

5 + e 2 

m v pa m I e r) I , (40) 






up to order one coefficients in the small entries. To obtain a large solar mixing angle, 
we need \r} — 8\ < e. In realistic models there is no reason for a cancellation between 
independent perturbations and thus we assume both 5 < e and 77 < e. 

Consider first the case 5 « e and rj < e. The solar mixing angle 8 12 is large but 
not maximal, as indicated by the MSW LA solution. We also have Am 2 atm pa 4m 2 , 
A ™Ln ~ Am at m e 2 and 

AmL • (41) 



If 1] « e and 5 e, we still have a large solar mixing angle and Am 2 sun w e 2 Am 2 



a4m ' 



as before. However m ee will be much smaller than the estimate in (j41j) . This is the 
case of the models based on the above mentioned U(1)f flavor symmetry that, at 
least in its simplest realization, tends to predict 5 ~ e 2 . In this class of models we 
find 



/ Am 



2 



below the sensitivity of the next generation of planned experiments. It is worth to 
mention that in both cases discussed above, we have 

\ Ami 



atm , 



which might be close to the present experimental limit if the oscillation frequency of 
the LA solution for solar neutrinos is in the upper part of its allowed range. 

If both 5 and rj are much smaller than e, the 12 block of m u has an approximate 
pseudo-Dirac structure and the angle 9 12 becomes maximal. This situation is typical 
of some models where leptons have U(l)p charges of both signs whereas the order 
parameters of U(l)p breaking have all charges of the same sign^Hl. We have two 
eigenvalues approximately given by ±m e. As an example, we consider the case 
where r] = and 5 « e 2 . We find sin 2 29 12 ~ 1 — e 2 /4, Am 2 un « m 2 e 3 and 




In order to recover the MSW LA solution we would need a relatively large value of e. 
But this is in general problematic because, on the one hand the presence of a large 
perturbation raises doubts about the consistency of the whole approach and, on the 
other hand, in existing models where all fermion sectors are related to each other, e 
is never larger than the Cabibbo angle. 

Summarising, within normal hierachical models there is enough flexibility to re- 
produce in a natural way the experimental frequencies and mixing angles. In partic- 
ular the lopsided matrix solution of the large atmospheric mixing, inspired by SU(5), 
where the large mixing arises from the charged lepton sector, can be extended rather 
naturally to also account for the solar sector and for the small 9 13 mixing angle. 



8.1. Semi- anarchy 



We have seen that anarchy is the absence of structure in the neutrino sector. Here 
we consider an attenuation of anarchy where the absence of structure is limited to 



the 23 sector. The typical texture is in this case: 



/ 5 e e 
e 1 1 



) 



(45) 



\e 1 1 



where 5 and e are small and by 1 we mean entries of 0(1) and also the 23 determinant 
is of 0(1). We see that this texture is similar to eq. (f3TJ|) when rj ~ 0(1). Clearly, 
in general we would expect two mass eigenvalues of order 1, in units of m, and one 
small, of order 5 or e 2 . This pattern does not fit the observed solar and atmospheric 
observed frequencies. However, given that the ratio r = (Am 2 un ) LA / Am 2 tm ~ 1/40 
is not too small, we can assume that the small value of r is generated accidentally, 
as for anarchy. We see that, if we proceed with the same change of basis as from 
eq. (}39|) to eq. (J47i|) . it is sufficient that by chance r\ ~ 6 + e 2 in order to obtain the 
correct value of r with large #23 and 612 and small #13 ~ e. The natural smallness of 
#13 is the main advantage over anarchy. We will come back to this class of models in 
a following section. 

9. Grand Unified Models of Fermion Masses 

We have seen that the smallness of neutrino masses interpreted via the see-saw 
mechanism directly leads to a scale A for L non- conservation which is remarkably 
close to Mqut- Thus neutrino masses and mixings should find a natural context 
in a GUT treatment of all fermion masses. The hierarchical pattern of quark and 
lepton masses, within a generation and across generations, requires some dynamical 
suppression mechanism that acts differently on the various particles. This hierarchy 
can be generated by a number of operators of different dimensions suppressed by 
inverse powers of the cut-off A c of the theory. In some realizations, the different powers 
of 1 / A c correspond to different orders in some symmetry breaking parameter v f arising 
from the spontaneous breaking of a flavour symmetry. In the next subsections we 
describe some simplest models based on SU(5) x U(1)f and on SO(10) which illustrate 
these possibilities'^. It is notoriously difficult to turn these models into fully realistic 
theories, due to well-known problems such as the doublet-triplet splitting, the proton 
lifetime, the gauge coupling unification beyond leading order and the wrong mass 
relations for charged fermions of the first two generations. Here we adopt the GUT 
framework simply as a convenient testing ground for different neutrino mass scenarios. 

9.1. Models Based on Horizontal Abelian Charges 

We discuss here some explicit examples of grand unified models in the framework 
of a unified SUSY SU(5) theory with an additional U(l)p flavour symmetry. The 
SU(5) generators act "vertically" inside one generation, while the U(1)f charges are 



different "horizontally" from one generation to the other. If, for a given interaction 
vertex, the U(l)p charges do not add to zero, the vertex is forbidden in the symmetric 
limit. But the symmetry is spontaneously broken by the VEV's Vf of a number of 
"flavon" fields with non- vanishing charge. Then a forbidden coupling is rescued but 
is suppressed by powers of the small parameters Vf/A c with the exponents larger for 
larger charge mismatch!^. We expect M GUT & Vf ^ A c ^ M Pi . Here we discuss 
some aspects of the description of fermion masses in this framework. 

In these models the known generations of quarks and leptons are contained in 
triplets and (i = 1,2,3) corresponding to the 3 generations, transforming 
as 10 and 5 of SU(5), respectively. Three more SU(5) singlets describe the RH 
neutrinos. In SUSY models we have two Higgs multiplets H u and Hd, which transform 
as 5 and 5 in the minimal model. The two Higgs multiplets may have the same or 
different charges. In all the models that we discuss the large atmospheric mixing 
angle is described by assigning equal flavour charge to muon and tau neutrinos and 
their weak SU(2) partners (all belonging to the 5 = (l,d c ) representation of SU(5)). 
Instead, the solar neutrino oscillations can be obtained with different, inequivalent 
charge assignments. There are many variants of these models: fermion charges can all 
be non-negative with only negatively charged flavons, or there can be fermion charges 
of different signs with either flavons of both charges or only flavons of one charge. 
We can have that only the top quark mass is allowed in the symmetric limit, or that 
also other third generation fermion masses are allowed. The Higgs charges can be 
equal, in particular both vanishing or can be different. We can arrange that all the 
structure is in charged fermion masses while neutrinos are anarchical. 



9.1.1. F(fermions)> 



Consider, for example, a simple model with all charges of matter fields being non- 
negative and containing one single flavon 9 of charge F= —1. For a maximum of 
simplicity we also assume that all the third generation masses are directly allowed 
in the symmetric limit. This is realized by taking vanishing charges for the Higgses 
and for the third generation components ^1°, ^3 and ^3. If we define F(\l/f ) = qf- 
(R = 10, 5, 1; i — 1, 2, 3), then the generic mass matrix m has the form 



m 



(ynXtf+tf yi2^ +q " yi3^ +q ? \ 

y 2X \<#+<£ y 22 \ q * +q *' y w X^ 
\ysM +q ? 2/32 A<^r y 33 X^'j 



(46) 



T 



where all the yij are dimensionless complex coefficients of order one and m u , = m 
mo and M arise by choosing (R, R') = (10, 10), (5, 10), (1, 5) and (1, 1), respectively. 
We have A = (0)/A c and the quantity v represents the appropriate VEV or mass 
parameter. The models with all non-negative charges and one single flavon have 
particularly simple factorization properties. For instance in the see-saw expression for 



m u = rn/^M x vtld the dependence on the q\ charges drops out and only that from qf 
remains. In addition, for the neutrino mixing matrix C/^, which is determined by m u 

I 5 5 1 

in the basis where the charged leptons are diagonal, one can prove that U{j ~ A Qi ~ q ^\ 
in terms of the differences of the 5 charges, when terms that are down by powers 
of the small parameter A are neglected. Similarly the CKM matrix elements are 
approximately determined by only the 10 charges ^3; V^ KM ~ A'^ 0- ^ '. If the 
symmetry breaking parameter A is numerically close to the Cabibbo angle, we can 
choose: 

(9 1 10 ,g 2 10 ) g3 1 °) = (3,2,0) , (47) 

thus reproducing V us ~ A, V c b ~ A 2 and V u b ~ A 3 . The same q}° charges also fix 
m u : m c : m t ~ A 6 : A 4 : 1. The experimental value of m u (the relevant mass values 
are those at the GUT scale: m = m^McuT) ^^) would rather prefer q\° = 4. Taking 
into account this indication and the presence of the unknown coefficients ~ 0(1) 
it is difficult to decide between q\° = 3 or 4 and both are acceptable. Of course the 
charges (g|°, q^ , q^ ) = (2, 1, 0) would represent an equally good choice, provided we 
appropriately rescale the expansion parameter A. Turning to the 5 charges, if we take 

(qlqlq!) = (b,o,o) b>0 , (48) 

together with eq. ()47|) we get the patterns : m s : ~ m e : : m T ~ X 3+b : A 2 : 
1. Moreover, the 22, 23, 32, 33 entries of the effective light neutrino mass matrix m v 
are all 0(1), thus accommodating the nearly maximal value of S23. The small non 
diagonal terms of the charged lepton mass matrix cannot change this. We obtain: 

(X 2b x b x b \ 2 

m v = \ X b 1 1 \ V -f (A,SA) , (49) 
\X b 1 1 / A 

where v u is the VEVs of the Higgs doublet giving mass to the up quarks and all the 
entries are specified up to order one coefficients. If we take v u ~ 250 GeV, the mass 
scale A of the heavy Majorana neutrinos turns out to be close to the unification scale, 
A ~ 10 15 GeV. 

If b vanishes, then the light neutrino mass matrix will be structure-less and we 
recover the anarchical (A) picture of neutrinos discussed in section 6.1. In a large 
sample of anarchical models, generated with random coefficients, the resulting neu- 
trino mass spectrum can exhibit either normal or inverse hierarchy. For down quarks 
and charged leptons we obtain a weakened hierarchy, essentially the square root than 
that of up quarks. 

If b is positive, then the light neutrino mass matrix will be structure-less only in 
the (2,3) sub-sector and we get semi-anarchical (SA) models, introduced in section 
7.2. In this case, the neutrino mass spectrum has normal hierarchy. However, unless 
the (2,3) sub-determinant is accidentally suppressed, atmospheric and solar oscillation 



frequencies are expected to be of the same order and, in addition, the preferred solar 
mixing angle is small. Nevertheless, such a suppression can occur in a fraction of semi- 
anarchical models generated with random, order one coefficients. The real advantage 
over the fully anarchical scheme is represented by the suppression in U e 3. 

Note that in all previous cases we could add a constant to qf, for example by 
taking (qf, gf, g|) = (2 + 6,2,2). This would only have the consequence to leave 
the top quark as the only unsuppressed mass and to decrease the resulting value of 
tan f3 = v u /vd down to A 2 mt/m&. A constant shift of the charges q\ might also provide 
a suppression of the leading v c mass eigenvalue, from A c down to the appropriate scale 
A. One can also consider models where the 5 and 5 Higgs charges are different, as in 
the "realistic" SU(5) model of ref. Also in these models the top mass could be 
the only one to be non- vanishing in the symmetric limit and the value of tan /3 can 
be adjusted. 



9.1.2. F(fermions) and F(flavons) of both signs 

Models with naturally large 23 splittings are obtained if we allow negative charges 
and, at the same time, either introduce flavons of opposite charges or stipulate that 
matrix elements with overall negative charge are put to zero. For example, we can 
assign to the fermion fields the set of F charges given by: 

fe 10 ^ 2 10 ^3 ) = (3,2,0) 

(qlqli) = (6,o,o) &>2a>o 

(QlQiQl) = (a,-a,0) . (50) 

We consider the Yukawa coupling allowed by U(l)F-neutral Higgs multiplets in the 5 
and 5 SU(5) representations and by a pair 9 and 9 of SU(5) singlets with F= 1 and 
F= —1, respectively. If b = 2 or 3, the up, down and charged lepton sectors are not 
essentially different than in the SA case. Also in this case the 0(1) off-diagonal entry 
of mi, typical of lopsided models, gives rise to a large LH mixing in the 23 block which 
corresponds to a large RH mixing in the d mass matrix. In the neutrino sector, after 
diagonalization of the charged lepton sector and after integrating out the heavy RH 
neutrinos we obtain the following neutrino mass matrix in the low-energy effective 
theory: 

(X 2b X b X b \ 9i2 

m v = \ X b 1 + A a A' a 1 + A a A /a 
\X b l + A a A ,a l + A a A /a / 

where A' is given by (0) / A c and A as before denotes the large mass scale associated 
to the RH neutrinos: A ^> v u ^. The 0(1) elements in the 23 block are produced 
by combining the large LH mixing induced by the charged lepton sector and the 
large LH mixing in mo- A crucial property of m u is that, as a result of the see-saw 
mechanism and of the specific U(1)f charge assignment, the determinant of the 23 



x (H), (51) 



block is automatically of 0(A a A' a ) (for this the presenc e of negative charge values, 
leading to the presence of both A and A' is essential EHES!! ) rp ne neu trino mass matrix 
of eq. (joTj) is a particular case of the more general pattern presented in eq. (J3*9]l , for 
5 « X 2b , e k ^ and r] m X a X' a . If we take A m A', it is easy to verify that the 
eigenvalues of m u satisfy the relations: 

m x : m 2 : m 3 = A 2(f, ~ a) : A 2a : 1 . (52) 

The atmospheric neutrino oscillations require m 2 ~ 10~ 3 eV 2 . The squared mass 
difference between the lightest states is of 0(A 4a ) m|, not far from the MSW solution 
to the solar neutrino problem if we choose a = 1. In general U e s is non- vanishing, of 
0(X b ). Finally, beyond the large mixing in the 23 sector, m v provides a mixing angle 
#12 ~ X b ~ 2a in the 12 sector. When b = 2a, as for instance in the case 6 = 2 and 
a = 1, the MSW LA solution can be reproduced and the resulting neutrino spectrum 
is hierarchical (H). 

Alternatively, an inversely hierarchical (IH) spectrum can be obtained by choosing: 

(gi 10 ,g 2 10 ,? 3 10 ) = (3,2,0) 

(qlqiqt) = (1,-1,-1) 
?s) = (-1,1,0) 

(q Hu ,q Hd ) = (0,1) • (53) 

Due to the non- vanishing charge of the H d Higgs doublet, in the charged lepton sector 
we recover the same pattern previously discussed. The light neutrino mass matrix is 
given by: 

/ A 2 1 1 \ 
m v = 1 A' 2 A' 2 (IH) . (54) 

V 1 A' 2 A' 2 J 

The ratio between the solar and atmospheric oscillation frequencies is not directly 
related to the sub-determinant of the block 23, in this case. 

A representative set of models is listed in table 2. The hierarchical and the 
inversely hierarchical models may come into several varieties depending on the number 
and the charge of the flavour symmetry breaking (FSB) parameters. Above we have 
considered the case of two (II) oppositely charged flavons with symmetry breaking 
parameters A and A'. It may be noticed that the presence of two multiplets 9 and 9 
with opposite F charges could hardly be reconciled, without adding extra structure 
to the model, with a large common VEV for these fields, due to possible analytic 
terms of the kind (99) n in the superpotential. Therefore it is instructive to explore 
the consequences of allowing only the negatively charged 9 field in the theory, case I. 
In case I, it is impossible to compensate negative F charges in the Yukawa couplings 
and the corresponding entries in the neutrino mass matrices vanish. Eventually these 
zeroes are filled by small contributions, arising, for instance, from the diagonalization 



Table 2: Models and their flavour charges. 



Model 


*10 


*5 




(H u , H d ) 


Anarchical (A) 


(3,2,0) 


(0,0,0) 


(0,0,0) 


(0,0) 


Semi-Anarchical (SA) 


(2,1,0) 


(1,0,0) 


(2,1,0) 


(0,0) 


Hierarchical (H) 


(3,2,0) 


(2,0,0) 


(1,-1,0) 


(0,0) 


Inversely Hierarchical (IH) 


(3,2,0) 


(Irlrl) 


(-i,+i,o) 


(o,+i) 



of the charged lepton sector or from the transformations needed to make the kinetic 
terms canonical. 

Another important ingredient is represented by the see-saw mechanism^. Hier- 
archical models and semi- anarchical models have similar charges in the (10, 5) sectors 
and, in the absence of the see-saw mechanism, they would give rise to similar results. 
Even when the results are expected to be independent from the charges of the RH 
neutrinos, as it is the case for the anarchical and semi-anarchical models, the see-saw 
mechanism can induce some sizeable effect in a statistical analysis. For this reason, 
for each type of model, but the hierarchical ones (the mechanism for the 23 sub- 
determinant suppression is in fact based on the see-saw mechanism), it is interesting 
to study the case where RH neutrinos are present and the see-saw contribution is 
the dominant one (SS) and the case where they are absent and the mass matrix is 
saturated by the non-renormalizable contribution (NOSS). 

With this classification in mind, we can distinguish the following type of models, 
all supported by specific choices of U(l) charges: Ass, Anoss, SAss, SAnoss, H(ss,i), 
H ( ss,ii), IH( SS ,i), IH( SS ,ii), IH(noss,i) and IH (NO ss,ii)- In table 3 we collect the predic- 
tions for the oscillation parameters arising from the neutrino mass matrices in eqs. 

(gHEEES}. 

It is interesting to quantify the ability of each model in reproducing the observed 
oscillation parameters. For anarchy, it has been observed that random generated, 
order-one entries of the neutrino mass matrices (in appropriate units), correctly fit 
the experimental data with a success rate of few percent. It is natural to extend this 



Table 3: Order of magnitude predictions for oscillation parameters, from neutrino mass matrices A, 
SA, H, IH in the text; c?23 denotes the sub-determinant in the 23 sector and we show the effect of 
its accidental suppression for the semi- anarchical model. In the estimates we have chosen A = A'. 
Inverse hierarchy predicts an almost maximal #12. 



Model 


parameters 


d 2 3 


Amj 2 /\Am 2 23 \ 




tan 2 812 


tan 2 623 


A 


b = 


0(1) 


0(1) 


0(1) 


0(1) 


0(1) 


SA 


b = 1 


0(1) 


0(4,) 


0(A) 


0(A 2 /d 2 3 ) 


0(1) 


Hn 


a = 1, b = 2 


0(A 2 ) 


0(A 4 ) 


0(A 2 ) 


0(1) 


0(1) 


Hi 


a = 1, b = 2 





0(A 6 ) 


0(A 2 ) 


0(1) 


0(1) 


IH 




0(A 4 ) 


0(A 2 ) 


0(A 2 ) 


1+0(A 2 ) 


0(1) 



analysis to include also the other models based on SU(5) x U(l)^*, which have mass 
matrix elements defined up to order-one dimensionless coefficients y^ (see eq. HBJ). For 
each model, successful points in parameter space are selected by asking that the four 
observable quantities 0\ — r = Am 2 2 /|Am2 3 |, O2 = tan 2 fl 12 , 3 = |[/ e 3| = | sin^i 3 | 
and O4 = tan 2 623 fall in the approximately 3a allowed ranges H ^U^ : 

0.01 < r < 0.2 

024 < tan 2 6 12 < 0.89 (LA) (55) 

0.33 < tan 2 623 < 3.3 

The coefficients of the neutrino sector are random complex numbers with absolute 
values and phases uniformly distributed in intervals X = [0.5, 2] and [0, 2ir] respec- 
tively. The dependence of the results on these choices can be estimated by varying X. 
For each model an optimization procedure selects the value of the flavour symmetry 
breaking parameter A = A' that maximizes the success rate. The success rates are dis- 
played in fig. [U and El separately for the NOSS and SS cases. The two sets of models 
have been individually normalized to give a total rate 100. Before normalization the 
total success rates for NOSS and for SS were in the ratio 1.7:1. The present data are 



most easily described by the IH schemes in their NOSS version. Their performances 
are better by a factor of 10-30 with respect to the last classified, the anarchical mod- 
els. The ability of the IH schemes in describing the data can be appreciated from 
the distributions of the four observables, which, for IH( NO ss,ii) an d A = A' = 0.25, are 
displayed in fig. El 




Figure 1: Relative success rates without see-saw. The sum of the rates has been normalized to 100. 
The results correspond to the default choice 2" = [0.5,2], and to the following values of A = A': 
0.2, 0.2, 0.25, 0.3 for the models AnosSj SAnoss> IH(noss,ii)i and IH( N oss,i)> respectively (in our 
notation there are no H(noss.i)j H(nossji) models). The error bars represent the linear sum of the 
systematic error due to the choice of X and the statistical error. 

The observables r and U e 3 are strongly correlated!^. Actually, as discussed in ref. 

and shown in table 2, in inversely hierarchical models U e 3 is typically of order r. 
Therefore, once A and A' have been tuned to fit r, this choice automatically provides 
a good fit to U e3 . Moreover tan 2 6*12 is peaked around 1. At present tan 2 # 12 = 1 
is excluded for the LA solution, but, thanks to the width of the distribution, the 
experimentally allowed window is sufficiently populated. The width of the distribution 
is almost entirely dominated by the effect coming from the diagonalization of the 
charged lepton sector. Indeed, by turning off the small parameters A and A' in the 
mass matrix for the charged leptons, we get a vanishing success rate for IH( N oss,i)> 
whereas the rate for IH(noss,ii) decreases by more than one order of magnitude. It 
is worth stressing that even a moderate further departure of the window away from 
tan 2 612 = 1 could drastically reduce the success rates of the IH schemes. Finally, the 
tan 2 623 distribution is rather flat, with a moderate peak in the currently favoured 
interval. All the IH models, with or without see-saw, have distributions similar to 
those shown in fig. El In particular the tan 2 623 distribution of fig. El is qualitatively 
common to all U(l) models. This reflects the fact that the large angle 823 is induced 
by the equal charges F(52) = F(5s), a feature shared by all the models we have 
considered. 
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Figure 2: Relative success rates with see-saw. The sum of the rates has been normalized to 100. The 
results correspond to the default choice X = [0.5, 2], and to the following values of A = A': 0.2, 0.3, 
0.35, 0.5, 0.15, 0.2 for the models A ss , SA SS , H( SS)II ), H(ss,i), IH(ss,n) and IH (SS j), respectively. The 
error bars represent the linear sum of the systematic error due to the choice of X and the statistical 
error. 

Without see-saw mechanism, the next successful model is the semi-anarchical 
model SA. Compared to the IH case, the tan 2 9n distribution has no pronounced 
peak. Possible shifts in the central value of tan 2 9 12 would not drastically modify 
the results for the SA model. The distribution is peaked around A = 0.2 with 
tails that exceed the present experimental bound. The r distribution is centered near 
r = 1. Finally, the anarchical scheme in its NOSS version is particularly disfavoured, 
due to its tendency to predict r close to 1 and also due to U e 3, that presents a broad 
distribution with a preferred value of about 0.5. 

The overall picture changes significantly if the light neutrino mass matrix arises 
from the see-saw mechanism, as illustrated in fig. |2j The IH models are still rather 
successful. Compared to the NOSS case, the IH SS models slightly prefer higher values 
of r and, due to a smaller A = A', they have very narrow tan 2 9u distributions. The 
other distributions are similar to those of figs. EJ Equally good or even better results 
are obtained by the H(ss,ii) model, with distributions shown in fig. HJ We see that, at 
variance with the IH models, tan 2 9 12 is not spiky, which results in a better stability 
of the model against variation of the experimental results. The preferred value of 
r is close to the lower end of the experimental window. The U e ^ distribution is 
nicely peaked around A 2 . The SAss model is significantly outdistanced from H(ss,ii), 
IH(ss,ii) and IH(ss,i). It is particularly penalized by the U e ^, distribution, centered 
around A = 0.3. Finally, the least favoured models are H(ss,i) and A ss . The model 



IH( NO ss,ii) 




Figure 3: Distributions for IH(noss,ii); 2" — [0-5, 2], A = A' = 0.25, obtained with 10000 points V. 



H(ss,i) fails both in U e $ (see fig. |3j), which tends to be too large for the preferred value 
of A = A' = 0.5 and in tan 2 # 12 . The Ass model, as its NOSS version, suffers especially 
from the U e % distribution (see fig. |SJ) which is roughly centered at 0.5, with only few 
percent of the attempts falling within the present experimental bound. A large U e $ 
can be regarded as a specific prediction of anarchy and any possible improvement 
of the bound on \U e ^\ will wear away the already limited success rate of the model. 
The distributions in tan 2 9 12 and tan 2 #23 are equally broad and peaked around 1. 
Compared to the NOSS case, Ass has a better r distribution, well located inside the 
allowed window. 

As a general comment we observe that our results are rather stable with respect 
to the choice of the interval I. With only one exception, namely the crossing between 
H(ss,i) and Ass, the relative position of the different models, according to their ability 
in describing the data, does not vary when we shrink X, or when we extend it to cover 



H(ssji) 




Figure 4: Distributions for H (SSJI) , 1 = [0.5, 2], A = A' = 0.35, obtained with 50000 points V. 




the full circle of radius 1 a . This stability would be partially upset by restricting 

a The results for X = [0, n] are independent on n, since changing n amounts to perform a renor- 
malization of all mass matrices by a common overall scale, which is not felt by the observables we 
have used in our analysis. 



to the case of real coefficients V. In that case the relative rates are comparable to 
those obtained in the complex case only for sufficiently wide intervals X, typically 
X « — VA] U [VA, l/VA]. If we further squeeze X around ±1 the rates of all 

models tend to zero. 

In conclusion, models based on SU(5) x U(l)p are clearly toy models that can 
only aim at a semiquantitative description of fermion masses. In fact only the order 
of magnitude of each matrix entry can be specified. However it is rather impressive 
that a reasonable description of fermion masses, now also including neutrino masses 
and mixings, can be obtained in this simple context, which is suggestive of a deeper 
relation between gauge and flavour quantum numbers. There are 12 mass eigenvalues 
and 6 mixing angles that are specified, modulo coefficients of order 1, in terms of 
a bunch of integer numbers (from half a dozen to a dozen), the charges, plus 1 
or more scale parameters. Moreover all possible type of mass hierarchies can be 
reproduced within this framework. In a statistically based comparison, the range 
of r and the small upper limit on U e s are sufficiently constraining to make anarchy 
neatly disfavoured with respect to models with built-in hierarchy. If only neutrinos 
are considered, one might counterargue that hierarchical models have at least one 
more parameter than anarchy, in our case the parameter A. However, if one looks 
at quarks and leptons together, as in the GUT models that we consider, then the 
same parameter that plays the role of an order parameter for the CKM matrix, for 
example, the Cabibbo angle, can be successfully used to reproduce also the hierarchy 
implied by the present neutrino data. 

9.2. GUT Models based on SO(10) 

Models based on SO(10) times a flavour symmetry are more difficult to construct 
because a whole generation is contained in the 16, so that, for example for U(1)f, one 
would have the same value of the charge for all quarks and leptons of each generation, 
which is too rigid. But the mechanism discussed so far, based on asymmetric mass 
matrices, can be embedded in an SO (10) grand-unified theory in a rather economic 

way ffi7teEl. The 33 entries of the U ion mass matrices can be obtained through 

the coupling I63I63IOH among the fermions in the third generation, I63, and a 
Higgs tenplet 10^. The two independent VEVs of the tenplet v u and Vd give mass, 
respectively, to t/v T and b/r. The key point to obtain an asymmetric texture is the 
introduction of an operator of the kind 16216^16316^ . This operator is thought to 
arise by integrating out an heavy 10 that couples both to I62I6// and to 16316^. 
If the 16h develops a VEV breaking SO(10) down to SU(5) at a large scale, then, 
in terms of SU(5) representations, we get an effective coupling of the kind 5 2 10 3 5#, 
with a coefficient that can be of order one. This coupling contributes to the 23 entry 
of the down quark mass matrix and to the 32 entry of the charged lepton mass matrix, 
realizing the desired asymmetry. To distinguish the lepton and quark sectors one can 



further introduce an operator of the form 16il6j 10,^45//, = 2, 3), with the VEV 
of the 45h pointing in the B — L direction. Additional operators, still of the type 
16il6jl6#16# can contribute to the matrix elements of the first generation. The 
mass matrices look like: 



m,. = 



(0 





m D = 
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e 



1 
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5' 
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+ e/3 
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-e/3 
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5' \ 


mi -- 
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a + e 


1 / 



Vd 



(56) 



(57) 



They provide a good fit of the available data in the quarks and the charged lepton 
sector in terms of 5 parameters (one of which is complex) . In the neutrino sector one 
obtains a large 9 23 mixing angle, sin 2 29i 2 ~ 6.6 • 1CT 3 eV 2 and 13 of the same order 
of 9i2- Mass squared differences are sensitive to the details of the Majorana mass 
matrix. 

Looking at models with three light neutrinos only, i.e. no sterile neutrinos, from 
a more general point of view, we stress that in the above models the atmospheric 
neutrino mixing is considered large, in the sense of being of order one in some zeroth 
order approximation. In other words it corresponds to off-diagonal matrix elements of 
the same order of the diagonal ones, although the mixing is not exactly maximal. The 
idea that all fermion mixings are small and induced by the observed smallness of the 
non diagonal Vqkm matrix elements is then abandoned. An alternative is to argue 
that perhaps what appears to be large is not that large after all. The typical small 
parameter that appears in the mass matrices is A ~ \jm~d[rn s ~ ^m^/rrir ~ 0.20 — 
0.25. This small parameter is not so small that it cannot become large due to some 
peculiar accidental enhancement: either a coefficient of order 3, or an exponent of the 
mass ratio which is less than 1/2 (due for example to a suitable charge assignment), or 
the addition in phase of an angle from the diagonalization of charged leptons and an 
angle from neutrino mixing. One may like this strategy of producing a large mixing 
by stretching small ones if, for example, he/she likes symmetric mass matrices, as 
from left-right symmetry at the GUT scale. In left-right symmetric models smallness 
of left mixings implies that also right-handed mixings are small, so that all mixings 
tend to be small, unless non-renormalizable mass operators with a suitable flavour 
pattern are introduced. Typically, by exploiting the freedom in the parameter space, 
in this set of models a large #23 may be accommodated. The large mixing angle for 
solar neutrinos requires however the introduction of ad hoc terms, like for instance 
higher- dimensional operators contributing to the light neutrino masses independently 



from the see-saw mechanism 



10. Conclusion 

By now there are rather convincing experimental indications for neutrino oscil- 
lations. The direct implication of these findings is that neutrino masses are not all 
vanishing. As a consequence, the phenomenology of neutrino masses and mixings is 
brought to the forefront. This is a very interesting subject in many respects. It is 
a window on the physics of GUTs in that the extreme smallness of neutrino masses 
can only be explained in a natural way if lepton number conservation is violated. If 
so, neutrino masses are inversely proportional to the large scale where lepton number 
is violated. Also, the pattern of neutrino masses and mixings interpreted in a GUT 
framework can provide new clues on the long standing problem of understanding the 
origin of the hierarchical structure of quark and lepton mass matrices. 

Neutrino oscillations only determine differences of mf values and the actual scale 
of neutrino masses remain to be experimentally fixed. The detection of Qv/3[3 decay 
would be extremely important for the determination of the overall scale of neutrino 
masses, the confirmation of their Majorana nature and the experimental clarification 
of the ordering of levels in the associated spectrum. The recent results from cosmol- 
ogy indicate that neutrino masses are not a major fraction of the cosmological mass 
density Q u ^ 1.5%. The decay of heavy right-handed neutrinos with lepton number 
non- conservation can provide a viable and attractive model of baryogenesis through 
leptogenesis. The measured oscillation frequencies and mixings are remarkably con- 
sistent with this attractive possibility. 

While the existence of oscillations appears to be on a ground of increasing so- 
lidity, many important experimental challenges remain. For atmospheric neutrino 
oscillations the completion of the K2K experiment, now stopped by the accident that 
has seriously damaged the Superkamiokande detector, is important for a terrestrial 
confirmation of the effect and for an independent measurement of the associated pa- 
rameters. In the near future the experimental study of atmospheric neutrinos will be 
further pursued with long baseline measurements by MINOS, OPERA, ICARUS. For 
solar neutrinos the continuation of SNO, KamLAND and the data from Borexino will 
lead to a more precise determination of the parameters of the LA solution. Finally a 
clarification by MINIBOONE of the issue of the LSND alleged signal is necessary, in 
order to know if 3 light neutrinos are sufficient or additional sterile neutrinos must be 
introduced, in spite of the apparent lack of independent evidence in the data for such 
sterile neutrinos and of the fact that attempts of constructing plausible and natural 
theoretical models have not led so far to compelling results. Further in the future 
there are projects for neutrino factories and/or superbeams aimed at precision mea- 
surements of the oscillation parameters and possibly the detection of CP violation 
effects in the neutrino sector. 



Pending the solution of the existing experimental ambiguities a variety of theoret- 
ical models of neutrino masses and mixings are still conceivable. Among 3-neutrino 
models we have described a number of possibilities based on degenerate, inverted hi- 
erarchy and normal hierarchy type of spectra. The normal hierarchy option appears 
to us as the most straightforward and flexible framework. In particular the large 
atmospheric mixing can arise from lopsided matrices. Then the observed frequencies 
and the large solar angle can also be obtained without fine tuning in models where 
the 23 subdeterminant is automatically suppressed. 

The fact that some neutrino mixing angles are large and even nearly maximal, 
while surprising at the start, was eventually found to be well compatible with a 
unified picture of quark and lepton masses within GUTs. The symmetry group at 
M GUT could be either (SUSY) SU(5) or SO(10) or a larger group. For example, we 
have seen that models based on anarchy, semianarchy, inverted hierarchy or normal 
hierarchy can all be naturally implemented by simple assignments of U(1)f horizontal 
charges in a semiquantitative unified description of all quark and lepton masses in 
SUSY SU(5)x U(l)p- Actually, in this context, if one adopts a statistical criterium, 
hierarchical models appear to be preferred over anarchy. 

All we know about neutrino masses is well in harmony with the idea and the mass 
scale of GUT's. As a consequence neutrino masses have added phenomenological 
support to this beautiful idea and to the models of physics beyond the Standard 
Model that are compatible with it. In particular, if we consider the main classes of 
new physics that are currently contemplated, like supersymmetry, technicolour, large 
extra dimensions at the TeV scale, little Higgs models etc, it is clear that the first one 
is the most directly related to GUT's. SUSY offers a well defined model computable 
up to the GUT scale and is actually supported by the quantitative success of coupling 
unification in SUSY GUT's. For the other examples quoted all contact with GUT's 
is lost or at least is much more remote. In this sense neutrino masses fit particularly 
well in the SUSY picture that sofar remains the standard way beyond the Standard 
Model. 
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